Speed of sound in a Bose-Einstein condensate 
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In the present work we determine the speed of sound in a Bose-Einstein condensate confined by 
an isotropic harmonic oscillator trap. The deduction of this physical parameter is done resorting 
to the TV-body Hamiltonian operator. The single-particle eigenfunctions that have been employed 
in this formalism are those stemming from the corresponding harmonic oscillator potential, and an 
expression for the dependence of this speed on the temperature is also deduced. These functions are 
used in the calculation of the scaterring length, etc. The situation for a Bose-Einstein condensate 
of sodium is evaluated and the corresponding speed of sound is obtained and compared against the 
known experimental outcomes. The possibility that the solution, to the existing discrepancy between 
experiment and theoretical predictions, could be given by the Zaremba-Nikuni-Griffin formalism is 
also explored. 
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I. INTRODUCTION 

A fundamental aspect in the context of quantum gases 
is related to the analysis of their collective excitations. 
Collective modes in uniform quantum fluids can be cate- 
gorized in two different realms, namely, collisionless and 
hydrodynamic modes. The former are associated to the 
dynamic self-consistent mean fields, whereas, the latter 
emerge as a consequence of the properties of the interac- 
tions The study of collective excitations could pro- 
vide a deep insight into some physical characteristics of 
these systems. In other words, the study of the exci- 
tations in a Bose-Einstein condensate (BEC) opens up 
a new realm in the comprehension of the structure of 
quantum fluids, at least in the regime of dilute quan- 
tum fluids. Clearly, in this direction there are several 
properties to be considered; among them we may find 
the role that trapping potentials and interactions among 
the constituents of the system play in the determination 
of certain features of a BEC. At this point let us focus 
our attention on the consequences of binary interactions 
upon some features of a BEC. In some cases, for instance, 
^"^ Kb or "^^Na, the size of the corresponding condensate 
suffers an enlargement due to the presence of the two- 
body repulsive forces. A further effect of these repulsive 
interactions is the fact that the central density of a non- 
ideal BEC at very low temperatures can be two or three 
orders of magnitude higher than the one associated to 
an ideal BEC 0, i.e., when the interaction among the 
particles of the gas is switched off. 

A knowledge of the speed of sound opens up several 
possibilities. For instance, the study of correlated mo- 
mentum excitations in the many-body condensate wave 
function [3, 01 • As a further comment let us mention that 
the dispersion relation for elementary excitations has a 
very different form when the possibility of an interaction 
among the particles is considered or not [S]. Since bi- 
nary collisions are not a frequent event in a BEC it is a 
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little bit surprising that the concept of interaction plays 
a primordial role in the determination of some physical 
features. The answer to this interrogant stems from the 
large coherent mean field associated to a BEC [6;] ■ These 
comments lead us to conclude that the comprehension of 
the speed of sound in a condensate defines an issue that 
bears an important physical relevance. In this context let 
us mention that it is believed that Gross-Pitaevskii mean 
field formalism captures the most essentialproperties as- 
sociated to the ground state of a BEC Q- Of course, 
the comparison of the theoretical predictions against the 
experimental outcomes provides a test of the validity of 
the mean field formalism, at least indirectly. 

Usually the deduction of the speed of sound in a BEC is 
done resorting to the linearized time-independent Gross- 
Pitaevskii equation in the so-called Tomas-Fermi limit. 
Afterwards, this expression is cast in the form of two 
quantum hydrodynamic equations (one for the density 
fiuctuations and the second one for the velocity). Fi- 
nally, the velocity is eliminated from the aforementioned 
equations and a differential equation for the density fluc- 
tuations is obtained [8] . Clearly, the aforementioned pro- 
cedure does not exhaust all the possible manners in which 
the speed of sound can be deduced. Another way starts 
from the N-body Hamiltonian in the second quantiza- 
tion formalism, introduces the Bogoliubov approximation 
and the resulting Hamiltonian is diagonalized by means 
of the Bogoliubov canonical transformation. The energy 
of the ground state of this diagonalized Hamiltonian al- 
lows us to calculate the pressure and speed of sound of 
the corresponding BEC [3]. A crucial point in all these 
approaches is the Mean Field Theory (MET) This 
formalism requires the use of several assumptions, one of 
them is related to the introduction of single-particle wave 
functions, a fact that can be tracked down to the Gibbs- 
Bogoliubov-Feynman equation [ll|, |l2|. The choice of 
these aforementioned single-particle eigenfunctions is a 
consequence of a minimization procedure, nevertheless, 
in some cases free-particles eigenfunctions Q, in other 
situations the eigenfunctions related to the particular 
trapping potential are employed 5], without proving if 
they are the minimizing case. 
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In the present work we determine the speed of sound 
in a Bose-Einstein condensate confined by an isotropic 
harmonic oscihator trap. The deduction of this physical 
parameter is done resorting to the iV-body Hamiltonian 
operator. The single-particle eigenfunctions that have 
been employed in this formalism are those stemming from 
the corresponding harmonic oscillator potential, and an 
expression for the dependence of this speed on the tem- 
perature is also deduced. These functions are used in 
the calculation of the scattering length, etc. The situ- 
ation for a Bose-Einstein condensate of sodium is eval- 
uated and the corresponding speed of sound is obtained 
and compared against the known experimental outcomes. 
Finally, we must mention that there is a discrepancy, 
the one emerges from the comparison between some of 
the extant theoretical predictions against the experimen- 
tal outcomes. Here we discuss also the possibility that 
the solution to this discrepancy could be given by the 
Zaremba-Nikuni-Grifhng formalism. 



II. MEAN FIELD THEORY AND SPEED OF 
SOUND IN A BEC 

Our starting point is the A^-body Hamiltonian opera- 
tor, in addition we assume that the collisions in the gas 
are, mainly, two-body interactions, this result is a con- 
sequence of the fact that we introduced, as additional 
condition, a dilute gas Under these restrictions the 
aforementioned Hamiltonian reads 
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In this last expression we have the following terms 



condensation are always, sufficiently, dilute. Let us ex- 
plain this fact; usually, particle separations, in the case 
of alkali atom vapours, have an order of magnitude of, 
approximately, lO^nm, whereas, the scattering length 
(here denoted by a) is two orders of magnitude smaller, 
a lOOflo, being oq the Bohr radius [Tj]. These com- 
ments entail that an {n + l)-body collision is less prob- 
able than an n-body collision, i.e., we may keep only 
two-body collisions. 

At this point no restriction upon the two-body interac- 
tion has been imposed. The next step, in this direction, 
concerns the introduction of some of the postulates of 
MET. Indeed, the assumption of very low temperature 
implies that the s-wave approximation can be used and 
the interatomic potential can be described by the so- 
called pseudo-potential. At low temperatures one of the 
features of two-body interactions is the emergence of the 
concept of scattering length as a fundamental idea 
This last comment does not allow us to evaluate ([3]). In- 
deed, we must know the set of single-particle functions, 
i.e., {u^(r)}. It has to be stressed that the introduc- 
tion of single-particle wavefunctions in ([T]) already im- 
plies the introduction, at least partially, of MET. This 
last assertion can be understood noting that MET as- 
sumes that the n-body symmetrized wavefunction can 
be replaced by a 1-body case where the bridge between 
these two cases is given by the deduction (via a mini- 
mization process stemming from the Gibbs-Bogoliubov- 
Eeynmann equation) of a set of single-particle wavefunc- 
tions which lead to the definition of the so-called Mean 
Eield Hamiltonian [Tl| . 

In the context of the calculations about the speed of 
sound in BEG usually two choices are made, namely, (i) 
The eigenfunctions of a three-dimensional harmonic os- 
cillator Q; (ii) Eree particle wavefunctions |9|]. These 
choices are made without checking if they correspond 
to the minimum required by the Gibbs-Bogoliubov- 
Eeynmann formalism. Of course, the first choice, for the 
case of a BEG trapped by a three-dimensional harmonic 
oscillator, seems to be a good conjecture since it reflects 
the symmetry of the trap. 
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III. SPEED OF SOUND AND NON- VANISHING 
TEMPERATURE 



<a,/3| F|7,e>= / ul{r}u}{r}V{r}u^{r}u,{r}d-'r[5) 



< a\ y(,)|/3 J <(r)V^(,)(f)^.^(r)dV. (4) 

Here V{r) denotes the two-body potential, whereas 
V(e) (r) depicts the trapping potential. Glearly, we have, 
explicitly, assumed that only two-body interactions are 
relevant for our case. An explanation for this approxima- 
tion is related to the fact that the systems employed in 



In the present work our trap will be an isotropic three- 
dimensional harmonic oscillator (with a frequency equal 
to id) and the corresponding eigenfunctions will be our 
choice for single-particle wavefunctions to be introduced 
in ([T|). Since we consider the limit of very low tempera- 
tures we assume that only the first excited state is popu- 
lated. In other words, ground state and first excited state 
are the only states populated. Eor an ideal BEG (no in- 
teractions among the particles of the gas) the number 
of particles in excited states (for temperatures below the 
condensation temperature) is a function of the trapping 

^ j , here the poten- 



potential, namely, = A^[l — 



3 



tial is represented by a [l3j. If we consider the presence 
of repulsive two-body interactions, then the number of 
particles in the ground and first excited state, iVo and 
A^e, respectively, for temperatures below the condensa- 
tion temperature (Tc), are given by 



3 8 W'N 



(5) 



where 7 is a real number. It is readily seen that this last 
expression does not vanish if the following two conditions 
are fulfilled: (i) 7 7^ 0, and (ii) a-1 = /3 + 1 ^a = 13 + 2. 
This last conclusion entails, since only the ground and 
first excited states are populated, that < a|a^a|/3 >= 
if a 7^ l3. These arguments allow us to deduce, very 
easily, that if jV'iv > denotes the wave function of our 
-/V-body system then 
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In this last expression a denotes the scattering length, 
and V the volume of the BEC. In the present case, in 
which there is no container of volume V, but a trap con- 
fines the gas, the definition of V requires a sound expla- 
nation. In this context we recall that a one-dimensional 
harmonic oscillator, whose frequency reads w, acting 
upon a quantum particle of mass m defines a length pa- 
rameter 
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The ground state wavefunction reads 
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Resorting to ([5]) we may cast this last expression in the 
following form 
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This last result provides us with an expression that 
shows a temperature dependence, for the case in which 
T < Tc- We now consider ([3]). In the calculation of this 
term we will assume that the involved particles are in the 
ground state. Clearly, energy has to be conserved. This 
condition entails that if we denote by a superindex the 
values of the occupation numbers before and after the 
collision, b and a, respectively, then energy conservation 
entails for our two involved particles 



The size of this one-dimensional system will be defined 
by the value x — xi such that 'ip{xi) = ip{x = 0)e~^, i.e., 
xi = V^l- 

Under these conditions a three-dimensional system has 
a volume given by ([IJ 
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V = 23/2^3. 



(9) 



With these arguments we now proceed to calculate ((IJ . 

The first term to be addressed involves the kinetic en- 
ergy. In order to do this we cast the operator as a 
function of the creation and annihilation operators [T^ 



In the case in which the eigenfunctions of single- 
particle are the free-particle functions, the evaluation of 
the corresponding integral is done resorting to momen- 
tum conservation ^ . In this direction we find our sought 
term 
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The corresponding term < a|V^|/3 > has to be cal- 
culated over all the possibilities for a and /3. Since, by 
hypothesis, we have assumed that only the ground and 
first excited states are populated, then < a|V^|/3 > does 
not vanish if and only if the following two conditions are 
fulfilled: (i) a or /3 represent the ground or first excited 
states; (ii) a = /3. This is no new restriction. Indeed, 
according to the rules that the annihilation and creation 
operators satisfy notice that < a|a''^a|/3 >= 7<5(q,+i^^_i), 
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Here r denotes the relative distance between the in- 
volved particles. This last expression leads us to the def- 
inition of the concept of scattering length 
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The usual situation is proportional to the integral of 
the interaction potential 13]. We now deduce the energy 
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of this system, as a function of the temperature, in the 
usual manner [9|] 
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The pressure {P{T)) reads 



3V Ftt 

4 



3- 



37V 



P{T) 
17 



-7V- 



2^2 



my5/3 3mF5/3 
27ra;i2 



TV 



Vtt 



T \ 3 

TV 



(16) 



a^iV 

'vV 



m 



y2 



3-2 



Vtt 



(17) 



The speed of sound {cs{T)) is given by 
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IV. DISCUSSION AND CONCLUSIONS 

Our last expression allows us to predict the speed of 
sound, as a function of the temperature, of course, this 
happens only in the regime T < Tc. Notice that our 
prediction does not include a dependence upon the am- 
plitude of the disturbance, a fact that matches with the 
experimental output 0, [IBl • 

It has already been recognized that many theoretical 
studies, in the realm of BEG, ignore the role that the 
thermally excited atoms play in the definition of the char- 
acteristics of the gas @. The assumption of vanishing 
temperature is not correct, as a matter of fact thermody- 
namics tells us that in a practical sense the achievement 
of this temperature is impossible [l^. Therefore, the 
deduction of the speed of sound under the assumption 
of T = is an approximation, the one should be im- 
proved. A point that has to be underlined in the present 
manuscript concerns this issue. Indeed, a fleeting glimpse 
at ([12]) tells us that the second term on the right hand-- 
side takes into account the contribution to the kinetic 
energy of the thermal cloud of the system as a function 
of the temperature, i.e., we do not assume T — 0. 

In the context of the assumptions here accepted, of 
course, we have, as mentioned before, discarded the pos- 
sibility of having, in this kinetic term, transitions be- 
tween excited states and the ground state. In order to 
have this case we must consider that not only the first 
excited state is populated, but also higher states. From 



this last comment we expect a very small contribution to 
the speed of sound stemming from this neglected possibil- 
ity. An additional simplification has been introduced, but 
now in relation with ^ . If we consider (jl4p we, imme- 
diately, notice that we have only considered interactions 
between particles in the ground state with particles in 
the same state. Of course, more possibilities are present, 
for example, a particle in the ground state might interact 
with a particle in the first excited state. Since the num- 
ber of particles in the ground state is much larger than 
those in excited states we expect that the probability of 
having a ground state-ground state interaction is larger 
than having a ground state-first excited state interac- 
tion, and this last one is larger than the first state-first 
state interaction. In this sense has to be understood the 
assumption. 

Let us now confront our theoretical prediction against 
the extant experimental results. We will consider the 
case of a BEG comprised by sodium atoms. This choice 
is done since we have already experimental results in this 
direction; iV = 5 x 10^ n = lO'^^m.-^, = 2 x lO'^K, 
m = 35.2X 10-27Kg, / ~ IQ-^m [H. Finally, the scatter- 
ing length, for sodium, has already been measured fl7| , 
namely, a = 2.75 x 10~^m. We need also an assumption 
for our temperature, here we assume T = Q.QTc- These 
values imply 



Cs = 2.2 X IQ-^m/s. 
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A careful look at the present measurement outputs en- 
tails that our result is not a bad one [l^, i.e., it provides 
the correct order of magnitude. The shortcomings of the 
present manuscript, at least in the realm of its compat- 
ibility with the experimental results, are also shared by 
other approaches [isl . [isj . The most intriguing interro- 
gant in this sense is related to the fact that in the region 
of the thermodynamical space in which the assumption of 
MFT is strongly satisfied the theoretical prediction has 
its worst behavior 18]. Several conjectures could be put 
forward, in order to solve this puzzle. The inclusion of 
additional terms, for instance, in connection with the dis- 
cussion about the number of significant terms related to 
([3]) would, surely, modify the result. In this direction an 
alternative way is related to the deduction of the speed of 
sound resorting to the generalized Gross-Pitaevski equa- 
tion (usually known as Zaremba-Nikuni-Griffin equation 
(ZNG) isl) the one takes into account a coupling be- 
tween the condensate and non-condensate components of 
the corresponding system. A physical motivation behind 
this statement can be found in the fact that sound can 
be understood as density waves. Glearly, changes in the 
density, which involve changes in the separation among 
particles, are determined by the interactions among the 
constituents of the gas. This explains in a very simple 
way why the idea of interaction plays a relevant role in 
the determination of the speed of sound. The ZNG equa- 
tion [g!] includes in the dynamics of a BEG a coupling 
between the condensate and the noncondensate part of 
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the gas. Wc may rephrase this last conjecture asserting 
that one of our assumptions, namely, no interaction be- 
tween condensate and noncondensate components could 
define a wrong premise in this context. In other words, 



the present approaches render the correct order of mag- 
nitude for this speed, but finer details involved in the 
correct deduction of it could stem from the use of the 
ZNG equation. 
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